New Measurement of the Scintillation Efficiency of Low-Energy Nuclear Recoils in 
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Particle detectors that use liquid xenon (LXe) as detection medium are among the leading tech- 
nologies in the search for dark matter weakly interacting massive particles (WIMPs). A key enabling 
element has been the low-energy detection threshold for recoiling nuclei produced by the interac- 
tion of WIMPs in LXe targets. In these detectors, the nuclear recoil energy scale is based on the 
LXe scintillation signal and thus requires knowledge of the relative scintillation efficiency of nuclear 
recoils, £ e ff- The uncertainty in C e s at low energies is the largest systematic uncertainty in the 
reported results from LXe WIMP searches at low masses. In the context of the XENON Dark 
Matter project, a new LXe scintillation detector has been designed and built specifically for the 
measurement of C e s at low energies, with an emphasis on maximizing the scintillation light detec- 
tion efficiency to obtain the lowest possible energy threshold. We report new measurements of C e s 
at low energies performed with this detector. Our results suggest a C c s which slowly decreases with 
decreasing energy, from 0.144 ± 0.009 at 15 keV down to 0.088±g;^g at 3 keV. 

PACS numbers: 95.35.+d, 14.80.Ly, 29.40.-n, 95.55.Vj 

Keywords: Liquid Xenon, Scintillation, Nuclear Recoil, Dark Matter 
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I. INTRODUCTION 



In recent years, liquid xenon (LXe) particle detectors 
[1-5] have achieved a large increase in target mass and 
a simultaneous reduction in backgrounds and are now 
among the leading technologies in the search for dark 
matter weakly interacting massive particles (WIMPs). 
The XENON100 experiment, with a target mass of 
62 kg and a measured electronic recoil background of 
< 10~ 2 events/kg/d/keV [6] is currently the most sensi- 
tive dark matter search in operation. Key to the perfor- 
mance of XENON100 and of future experiments based 
on the same principle, is the ability to detect low-energy 
recoiling nuclei in LXe. Since WIMPs are expected to 
interact primarily with atomic nuclei, the nuclear recoil 
energy scale is based on the LXe direct scintillation sig- 
nal and thus requires knowledge of the scintillation yield 
of nuclear recoils. The uncertainty in the nuclear re- 
coil energy scale at low energies is the largest system- 
atic uncertainty in the reported results from LXe WIMP 
searches [2, 7]. 

The scintillation yield, defined as the number of pho- 
tons produced per unit energy, depends on the type of 
particle as well as the energy [8, 9]. Since a precise mea- 
surement of the absolute scintillation yield is rather diffi- 
cult, the relative scintillation efficiency of nuclear recoils, 
jCcffj is the quantity that is used to convert the scintil- 
lation signals of LXe dark matter detectors in nuclear 
recoil energies. £ e ff is defined as the ratio of the scintil- 
lation yield of nuclear recoils to that of electronic recoils 
from photoabsorbed 122 keV 7 rays from a 57 Co source, 



at zero electric field. The scintillation yield of electronic 
recoils is also an energy-dependent quantity [10], but the 
122 keV 7 rays provide a reference with which scintilla- 
tion yields at other energies or of other particles can be 
compared. 

Two methods have been used to measure £ ff at differ- 
ent energies, an indirect method with the full spectrum 
comparison of the simulated response to the measured 
response from the irradiation with a neutron source, and 
a direct method with monoenergetic neutron fixed-angle 
scatters. 

Indirect measurements [11, 12] infer the energy depen- 
dence of C e g by comparing experimental data obtained 
with a neutron source and a Monte Carlo simulation of 
the expected nuclear recoil energy spectrum. In this way, 
any neglected factors will be absorbed in the energy de- 
pendence of £ G ff. Such factors can include uncertain- 
ties in the energy spectrum of the neutron source, ef- 
ficiency losses near threshold, energy dependence of se- 
lection cuts, etc, and are typically difficult to measure 
precisely. 

Direct measurements [13-19] are performed by record- 
ing fixed- angle clastic scatters of monoenergetic neutrons 
tagged by organic liquid scintillator detectors. The recoil 
energy of the Xe nucleus is then entirely fixed by kine- 
matics and approximately given by 



m n M Xc 
In + M Xc 
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where E n is the energy of the incoming neutron, m n and 
Mxe are the masses of the neutron and Xe nucleus, re- 
spectively, and 9 is the scattering angle. The spread in 
measured recoil energies mostly comes from the energy 
spread of the neutron source and the angular acceptance 
of the LXe and neutron detectors due to their finite sizes. 
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Since the start of the XENON dark matter project, our 
group has already performed two direct measurements of 
jC e g. In this article we report results from a new measure- 
ment of this quantity, with an improved apparatus and 
elaborated control of systematic uncertainties. These are 
the most accurate measurements of £ e ff in LXe to-datc, 
down to 3 keV recoil energy. 



II. EXPERIMENTAL SETUP 

The measurement of £ e ff was carried out by irradiat- 
ing a LXe detector with approximately monocnergetic 
neutrons produced by a sealed-tube neutron generator 
[26] and placing two organic liquid scintillator detectors 
at different azimuthally symmetric positions with respect 
to the neutron generator-LXe detector axis, see Fig. 1. 




FIG. 1: Schematic of the experimental setup. A sealed-tube 
neutron generator, where deuterons of energy Ed are incident 
upon a titanium deuteride target, produces neutrons at vari- 
ous angles <p. Some of the neutrons emitted at an angle tp = -| 
scatter in the LXe detector at an angle 9 and are tagged by 
two EJ301 organic liquid scintillators neutron detectors. 



lation efficiency of nuclear recoils at zero held. A poly- 
tetrafluoroethylene (PTFE) frame serves as a mounting 
structure and alignment guide for the PMTs so that each 
PMT window covers a side of the cubic active volume. 
The PMT assembly is housed in a stainless steel detector 
vessel, surrounded by a vacuum cryostat. The detector 
vessel has a special cross shape that closely follows the 
contours of the PMT assembly to minimize the proba- 
bility that neutrons scatter before or after an interaction 
in the active volume. The PTFE mounting structure is 
suspended from the top by a stainless steel rod fixed to 
a linear displacement motion feedthrough. The motion 
feedthrough allows the adjustment of the vertical posi- 
tion of the assembly from the outside. A schematic of 
the detector is shown in Fig. 2. 
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FIG. 2: Explosion drawing of the LXe detector. 



The emphasis for the design of the LXe detector has 
been placed on the reduction of non-active LXe and 
other materials in the immediate vicinity of the active 
LXe volume, and on the maximization of the scintil- 
lation light detection efficiency The active LXe vol- 
ume is a cube with sides of length 2.6 cm viewed by 
six 2.5 cm x 2.5 cm Hamamatsu R8520-406 Sel photo- 
multiplier tubes (PMTs), the same type as used in the 
XENON100 experiment [2] but selected for high quan- 
tum efficiency (QE). The PMTs have a special Bialkali 
photocathode for low temperature operation down to 
— 110° C and have an average room temperature QE of 
32% [27] at 178 nm, the wavelength at which Xe scintil- 
lates. The PMTs are operated with positive high voltage 
bias so that the PMT metal body and photocathode are 
at ground potential. This ensures that the active volume 
remains at zero electric field. This is a prerequisite to this 
measurement, since £ e ff is defined as the relative scintil- 



A schematic of the gas handling, liquefaction and re- 
circulation system developed for this experiment is shown 
in Fig. 3. The Xe is purified in the gas phase by circu- 
lating it through a hot getter and re-liquified efficiently 
using a heat exchanger [20] . The LXe temperature is kept 
constant with an Iwatani PDC08 pulse tube refrigerator 
(PTR) delivering 24 W of cooling power at 165 K with 
an air-cooled 1.5 kW helium compressor. Thanks to the 
PTR stability and low maintenance, the cooling system 
developed for this detector has enabled us to acquire data 
over a two month long, uninterrupted run. During the 
measurements, the LXe temperature was maintained at 
—94° C which corresponds to a vapor pressure of 2 atm. 
Further details on the design and performance of the cool- 
ing system developed for this experiment are presented 
in Ref. [20]. 

Neutrons with an average energy of 2.5 McV were pro- 
duced via the 2 H(d, n) 3 He reaction in a compact sealed- 
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tube neutron generator provided by Schlumbcrgcr. The 
generator was operated at deuteron energies of 60, 65, 
75, or 80 kcV and with deuterium beam currents rang- 
ing from 60 to 100 /iA. The target is a self-regenerating 
titanium deuteride (TiD) thick target. 



FIG. 3: Xenon gas system used for continuous purification of For non-relativistic deuterons (E d < 20MeV), the en- 

the LXe. ergy of neutrons emitted in the 2 H(d, n) 3 He reaction is 

given by [21]: 



pi/2 _ (m d m n E d ) 1/2 {m d m n E d cos 2 tp + (m He + %) [m He Q + (™-Hc - m d ) E d }} 1/2 

il/^y — - cos (p ~\~ \^2i j 



m Ho + m n m Hc 



where m d and toho are the deuteron and 3 He nucleus 
masses, respectively, Q is the Q-vahie of the reaction, 
and ip is the neutron emission angle. At small and large 
emission angles the neutron energy depends significantly 
on the deuteron energy. However, there is a minimum 
in both dE n /d<p and dE n /dE d at ip ~ 100°, and conse- 
quently the energy spread of neutrons produced is min- 
imal near this angle. For this reason, the neutron gen- 
erator was operated in a configuration where deuterons 
are accelerated vertically and where the neutrons inci- 
dent on the LXe detector are those produced at <p = ? 
(see Fig. 1). 

The expected neutron yield from the generator can be 
computed from its operating conditions, specifically the 
deuteron energy, E d , the deuterium beam current, and 
the monoatomic/diatomic deuterium beam composition. 
Since deuterons are stopped in the target, the neutron 
angular yield Y tot (E d ,ip) is calculated from 

Y tot (E d , ip) = xY(E d , ip) + 2yY(E d /2, ip) (3) 



with 



Y(E d ,<p) 



<\>n d 



-i -l 



dE' d 



(4) 



and the neutron energy distribution at angle tp from 



N tot (E n , ip) dE n = [xN(E d , ip) dE d 



+2yN(E d /2, ip) dE d ] 



dE n 
dE d 



(5) 



with 



N(E dl ip) dE d 



>n d 



o(E d ,ip) 



dE 
d.r 



(6) 

where <j> is the incident deuteron flux, n d the num- 
ber density of deuterium atoms in the target and p its 
density, a(E d ,ip) is the 2 H(d, n) 3 He differential cross- 
section, (dE/dx) d the deuteron stopping power of the 
target, and x/y the monoatomic/diatomic beam frac- 
tions. The 2 H(d, n) 3 He neutron production cross-section 
is taken from the ENDF/B-VII.O database [22], the stop- 
ping power for protons from the PSTAR database [28], 
and the monoatomic/diatomic deuterium beam fraction 
is taken as 0.05/0.95 [29]. The result of the computation 
is shown in Fig. 4. 

Measurements of the neutron flux were also carried 
out using a Nuclear Research Corporation NP-2 portable 
neutron monitor. The NP-2 neutron monitor is a boron 
trifluoride proportional counter surrounded by a low den- 
sity polyethylene cylinder to moderate the neutron flux 
before capture by boron. Its accuracy is ±10%. The re- 
sult of the computation show remarkable agreement with 
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the measurement as a function of applied high voltage 
(maximum deuteron energy). The agreement with the 
measurement as a function of beam current is also good 
but worsens as the beam current increases. This effect is 
expected since, as the measured beam current increases, 
a larger fraction of it can be due to electronic leakage cur- 
rent instead of deuterium current. The energy spread of 
neutrons is assumed to be dominated by the energy loss 
of deuterons in the target and by neutrons produced at 
different angles that scatter back into the LXe detector 
direction. 
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FIG. 4: Measurement (points) and theoretical calculation 
(solid line) of the neutron generator flux as a function of high 
voltage (top) and beam current (down). 



The organic liquid scintillator neutron detectors used 
in this measurement are 3 inch diameter Eljen Tech- 
nologies M510 detector assemblies filled with EJ301 [30]. 
The liquid scintillator is encapsulated in an aluminium 
cylindrical container, coupled to a 3 inch ET Enterprises 
982IB PMT. EJ30I is especially adapted to fast neutron 
detection in the presence of 7 radiation due to its excel- 
lent pulse shape discrimination (PSD) characteristics. 

The alignment of the LXe detector and the neutron 
generator was performed with an auto-levelling laser 
mounted on a tripod. The horizontal line and the verti- 
cal line projected by the laser were used to set the height 
of the neutron generator target to the center of the LXe 
detector and to mark the position of both on the labora- 
tory floor. The laser was also used to align the two EJ30I 
neutron detectors with respect to the LXe detector. For 
each angle, the distance between the LXe detector and 
the two EJ301 neutron detectors was chosen to produce a 
recoil energy spectrum with a spread due to the angular 



acceptance of the cell of 10% to 20%. The desired EJ301 
neutron detector positions were marked on the floor us- 
ing 1.5 m aluminium rules, while their height was set 
with the help of the laser and a vertical rule. The EJ301 
neutron detectors were supported at the desired positions 
around the LXe detector by their own laboratory stands. 

The horizontal distance between the neutron generator 
and the LXe detector was fixed at 40 cm. The distance 
from the LXe detector to the EJ301 liquid scintillators 
varied from 100 cm, for the scattering angles correspond- 
ing to low-energy recoils, to 40 cm, for the scattering an- 
gles corresponding to higher energies. The positioning 
accuracy of the EJ301 neutron detectors is estimated to 
be better than 5 mm. 

The data acquisition electronics are largely the same 
as those of the XENON100 experiment. The signals from 
the six LXe PMTs arc fed into a Phillips 776 xlO am- 
plifier with two amplified outputs per channel. The first 
output of each channel is digitized by a 14-bit CAEN 
V1724 100 MS/s flash ADC with 40 MHz bandwidth 
while the second output is fed to a Phillips 706 lead- 
ing edge discriminator. The discriminator thresholds are 
set at a level of -15 mV, which corresponds to 0.5 photo- 
electrons (pe). The logic signals of the six discriminator 
outputs are added with a linear fan-in and discriminated 
to obtain a 2-fold PMT coincidence condition. The 2-fold 
PMT coincidence logic signal is passed to a 10 [is hold- 
off circuit to prevent re-triggering on the tail of the LXe 
scintillation signal and constitutes the LXe trigger. 

The signals from the two EJ301 PMTs are also fed into 
the x 10 amplifier, with one output digitized by the flash 
ADC unit, and the other output discriminated. However, 
since we do not expect a signal for both EJ301 PMTs in 
the same event, the two signals are multiplexed into a 
single channel before digitization. The two discriminator 
outputs are delayed and also digitized to serve as a de- 
multiplexing code. A copy of the discriminator outputs is 
combined into an OR gate and forms the EJ301 trigger. 

The neutron time-of-flight (TOF) is measured with an 
Ortec 566 time to amplitude converter (TAC), where the 
LXe trigger is used as the "start" signal and a delayed 
copy of the EJ301 trigger as the "stop" signal. The delay 
of the copy of the EJ301 trigger signal can be varied to 
calibrate the TOF measurement. The TAC output signal 
is also digitized by the flash ADC unit and the TOF is 
computed by the event processing program. 

Finally, for neutron scattering measurements, the trig- 
ger is taken as a coincidence within a 200 ns window of 
the LXe trigger and the EJ301 trigger. The low-energy 
roll-off of the trigger efficiency can contribute signifi- 
cantly [23] to the systematic uncertainty on £ ff, if it is 
not properly understood. As an improvement over previ- 
ous measurements of £ e ff , where the efficiency was based 
solely on simulations, the efficiency of the trigger setup 
described above was measured in addition to its simula- 
tion. The efficiency measurement was done using a 22 Na 
source (a /3 + emitter) placed between the LXe detector 
and a sodium iodide Nal(Tl) detector. The back-to-back 
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pair of 511 keV 7 rays from the f3 + annihilation interact 
effectively at the same time in the LXe and the Nal(Tl) 
detectors. The Nal(Tl) detector was positioned such that 
the solid angle it subtended at the source was larger than 
the one subtended by the active LXe volume. This en- 
sures that the whole active volume of the LXe detector 
is probed. For this measurement, the triggering signal 
consisted of the discriminated signal of the output of the 
Nal(Tl) detector. In addition to the signals of the LXe 
PMTs, the trigger signal of the normal LXe trigger was 
digitized with the flash ADC. The efficiency is inferred 
by computing the fraction of events accompanied by a 
LXe trigger signal as a function of their measured num- 
ber of photoelectrons. Fig. 5 shows the resulting trigger 
efficiency 

The expected efficiency has also been computed via a 
detailed Monte Carlo simulation. The simulation takes 
into account the spatial distribution of the light collec- 
tion efficiency within the active LXe volume, PMT quan- 
tum efficiencies, PMT first dynode collection efficiencies, 
PMT gains, variations in the single photoelectron pulse 
height distribution, and discriminator trigger thresholds. 
The result is also shown in Fig. 5. The discrepancy 
between the simulated and measured efficiencies is at- 
tributed to the limited accuracy in modelling hardware 
components of the trigger. Consequently, the measured 
trigger efficiency is used to extract £ e ff (Sec. Ill E) from 
the neutron scattering measurements. In effect, this dis- 
crepancy only stresses the need to perform a measure- 
ment of the trigger efficiency instead of relying solely on 
simulation. 
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FIG. 5: Measured (points) and simulated (solid curve) effi- 
ciency of the 2-fold coincidence LXe trigger. The measured 
trigger efficiency is used to extract £ e fl (Sec. HIE) from the 
neutron scattering measurements. 



III. DATA ANALYSIS 
A. Calibration 

The LXe PMT gains were measured under single pho- 
toelectron conditions with a pulsed blue light emitting 
diode (LED), embedded in the PTFE mounting struc- 
ture. The gains were equalized to a mean value of 
2.0 x 10 6 by adjusting the individual PMT anode bias 
voltages at the beginning of the experiment and were 
monitored regularly throughout. At the operating con- 
ditions chosen, the single photoelectron resolution of the 
LXe PMTs has a mean value of 60%. 

The electronic recoil energy scale is calibrated using 
122 keV 7 rays since £ e ff is defined as the scintilla- 
tion light yield of nuclear recoils relative to the yield 
of 7 rays of that energy. Calibrations with an exter- 
nal 100 /iCi 57 Co source were taken regularly through- 
out the duration of the experiment. Fig. 6 shows the 
scintillation spectrum of a calibration with the 57 Co 
source. The scintillation light yield was measured to be 
L y = 24.14±0.09(stat)±0.44(sys)pe/keVee with a reso- 
lution (a/E) of 5%. This very high light yield, combined 
with the 90% trigger efficiency at 7 pe, implies that en- 
ergy spectra do not suffer from efficiency losses down to 
energies as low as 0.3 keVee. 

Since the attenuation length of 122 keV 7 rays in 
LXe is 3 mm [24] the external 57 Co source mostly probes 
the light yield in the outer layers of the active volume. 
However, since the elastic scattering mean free path of 
2.5 MeV neutrons in LXe is ~ 20 cm, the expected spa- 
tial distribution of nuclear recoils is uniform. Hence, the 
quantity of interest is the average light yield over the 
whole volume. The spatial uniformity of the light col- 
lection has been investigated using a light propagation 
simulation that takes into account the detailed geometry 
of the active volume, the photocathode coverage of the 
PMTs, and the reflectivity of PTFE. The relative light 
yield variation over the active volume was calculated to 
be less than 2%, with a maximum variation of 5% near 
the edge. 

The TOF calibration was performed with the 22 Na 
source placed between the LXe detector and the EJ301 
neutron detectors. The two emitted 511 keV 7 rays emit- 
ted interact in both detectors and give a TOF = cal- 
ibration point. The time delay between the LXe and 
EJ301 triggers was varied with a delay generator over a 
range of 32 ns to obtain the TOF calibration. 

B. Data Processing, Selection 

Events consist of the six LXe PMT waveforms, the 
waveform of the triggered EJ301 neutron detector, as well 
as the waveform of the TAC output. For each event, the 
processing software searches the LXe PMT waveforms 
for scintillation signals and records several parameters 
for each pulse found: position, area, height, etc. The 
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FIG. 6: Scintillation light spectrum of the 100 /xCi 57 Co 
source used to calibrate the LXe light yield, L v . The 
peak at ~ 3000 pe is the 122 keV photoelectric absorption 
peak. Also visible is the Xe 30keV characteristic X-ray at 
~ 700 pe. This calibration gives a scintillation light yield of 
L y = 24.3pe/keV. 



FIG. 7: Neutron interactions (gray/black) and gamma inter- 
actions (violet) for one of the two EJ301 neutron detectors. 
The neutron events selected by the PSD cut (red dashed line) 
and the EJ301 energy threshold cut (blue dashed line) are 
shown in black. 



same procedure is applied to the EJ301 neutron detector 
waveform. In addition, a PSD parameter is defined as the 
fraction of the total EJ301 scintillation signal contained 
within the tail of the pulse. The tail of the pulse is defined 
as a region with a lower boundary at 30 ns after the pulse 
peak and an upper boundary at the time where the pulse 
returns to 1% of the peak amplitude. Finally, the TOF 
is computed from the height of the TAC ouput and the 
TOF calibration values. For a LXe detector and EJ301 
neutron detector distance of 1 m, the the typical neutron 
TOF is 45 ns. 

Three selection cuts are applied to the data: a PSD 
cut to select neutron interactions in the EJ301 neutron 
detectors, a lower energy threshold cut on the EJ301 scin- 
tillation signal, and a TOF cut. 

The PSD cut uses the j/n discrimination capabilities 
of the EJ301 liquid scintillator and allows the selection 
of neutron interactions in the EJ301 with high efficiency. 
However, since this cut is based on the scintillation sig- 
nal, the discrimination power degrades with decreasing 
energy and it is preferable to apply an energy threshold 
cut on the EJ301 scintillation signal. This reduces the 
background from neutrons that scattered in other ma- 
terials, in addition to the scatter in the LXe detector, 
since the neutron energy was reduced and hence cannot 
deposit the maximum energy expected. Fig. 7 shows the 
events selected by the combination of the PSD cut and 
the EJ301 energy threshold cut for one of the two EJ301 
neutron detectors. 

A subtle effect can be observed in the TOF calibration 
data. Since the LXe trigger comes from the coincidence 
of two LXe scintillation photons, the LXe scintillation 
light decay time systematically shifts the TOF measure- 



ment of scintillation signals with fewer photons to lower 
TOF values. The effect is more pronounced for 7 rays 
interacting in LXe due to the slower recombination time 
for electron recoils. The TOF value is corrected for this 
effect in the event processing software. However, the cor- 
rection does not eliminate the spread in TOF values due 
to this effect. To eliminate any possible bias, we chose 
TOF cuts that contain the full neutron TOF peak for 
each scattering angle measurement, even if this results 
in a higher contamination of the recoil spectrum by neu- 
trons that scattered in other materials. 



C. Measured Recoil Distributions 

Neutron scattering data were acquired at eight differ- 
ent angles: 23°, 26.5°, 30°, 34.5°, 39.5°, 45°, 53°, and 
120°, corresponding respectively to recoil energies E nl of 
3.0±0.6, 3.9±0.7, 5.0±0.8, 6.5±1.0, 8.4±1.3, 10.7±1.6, 
14. 8 ± 1.3, and 55.2 ±8.8 keV. The uncertainty in nuclear 
recoil energies, dominated by the angular acceptance of 
the detectors, is extracted from the results of the Geant4 
Monte Carlo simulation described in Sec. HID. For most 
angles, the first EJ301 neutron detector was placed in 
the plane of the neutron generator and the LXe detector, 
while the second was placed at the same distance to the 
LXe detector but higher above the laboratory floor, with 
an azimuthal angle of about 45°. For all measurements 
the recoil energy and TOF distributions of the two EJ301 
neutron detectors are compatible. Figures 8 and 9 show 
the measured recoil energy and TOF spectra for the all 
scattering angles. 

All measurements show two very clear peaks in the 
TOF spectrum. The peak at TOF = corresponds to 7 
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rays that Compton scatter in the LXe detector before in- 
teracting in the EJ301 neutron detectors, while the peak 
at later TOF values corresponds to neutrons. For the 
lower recoil energies the effect of the spread of the TOF 
distribution to smaller values, discussed in the previous 
section, is noticeable in Figures 8 and 9. 

For energies of 6.5 keV and above, the peak in the 
recoil spectrum is clearly above the beginning of the low- 
energy trigger efficiency roll-off. For these energies, £ e g 
could even be computed directly without much uncer- 
tainty, simply by fitting a Gaussian to the peak. For 
energies below 6.5 keV, a more sophisticated procedure 
that takes into account the trigger efficiency is warranted. 
The procedure used to extract C e s from the measured re- 
coil energy spectra is detailed in Sec. HID. 



D. Monte Carlo Simulation 

Extensive Geant4 simulations of the expected neutron 
scattering rate, nuclear recoil energy distribution and 
neutron TOF distribution were performed for each scat- 
tering angle. Each simulation takes into account a real- 
istic description of the neutron generator, LXe detector, 
detector vessel, vacuum cryostat, support frame, as well 
as the measured positions of the two EJ301 neutron de- 
tectors. The livetime of each scattering angle simulation 
is calculated from the expected neutron yield of the gen- 
erator at the operating conditions of the measurement 
and the results are scaled accordingly. 

The information recorded in the simulation includes 
the energy, position, time, type of particle, and physical 
process responsible for each energy deposit in the LXe de- 
tector, as well as the energy, time, and type of particle for 
each energy deposit in the EJ301 neutron detectors. The 
energy and direction of the primary neutron is sampled 
from the calculated energy-angle distribution of neutrons 
produced by the generator, as given by eq. 6. 

Neutrons that interact in the active LXe volume can 
deposit energy via elastic or inelastic scattering, once or 
multiple times, and may additionally scatter in materials 
outside of the active volume. The contributions from all 
these classes of events has been inferred from the Monte 
Carlo simulation. Figs. 10 and 11 show the simulated re- 
coil energy and TOF distributions for all measured scat- 
tering angles. The energy spectrum of elastic recoils con- 
sists of a peak, roughly centered at the recoil energy cor- 
responding to the angle at which the EJ301 neutron de- 
tectors have been placed, and an approximately exponen- 
tially distributed background. The peak is due to neu- 
trons that clastically scatter once in the LXe detector and 
interact nowhere else (pure single elastic scatters), while 
the exponential background is due to neutrons that addi- 
tionally scatter in other materials surrounding the active 
volume, that is, their energy deposit in the LXe is essen- 
tially that of a recoil with a random scattering angle. For 
the smaller angles the nuclear recoil energy spectrum is 
clearly dominated by pure single elastic recoils with 61% 



at 23° (3 keV) but the proportion gradually decreases 
reaching 49% at 8.4 keV and 17.8% at 55.2 keV. The 
multiple scatter contribution is negligible at all angles, 
from 3.7% at 23° up to 7.0% at 120°, due to the small 
dimensions of the active LXe volume compared to the 
neutron elastic scattering mean free path. These simu- 
lation results demonstrate clearly that the design goal of 
minimizing the amount of materials in the vicinity of the 
active LXe volume has been achieved. 



E. Extracting £ e fi 

The energy dependence of £ e ff is obtained by minimiz- 
ing the x 2 statistic between the measured recoil distribu- 
tion and the simulated distribution with respect to two 
free parameters, C e fi,j = £cfi(E nr j) and Rj = R(E nr j), 
respectively the scintillation efficiency and the energy res- 
olution at the recoil energy measured E nr j. These fits 
are shown in Fig. 12 for all measured scattering angles. 
Explicitly, the x 2 statistic is computed from 

where hi and gi are the measured and simulated event 
rates in energy bin i, respectively, and Oh,i and a g ^ the 
uncertainties in the measured and simulated event rates 
in bin i, respectively. The bins over which the \ 2 statistic 
is computed varies depending on the scattering angle so 
that the x 2 does not become dominated by effects in the 
higher energy tail of the recoil distribution, nor the low 
trigger efficiency region. Bins with no counts are removed 
from the sum (and the number of degrees of freedom is 
decreased by one). 

The steps involved in transforming the simulated re- 
coil energy distributions of Figs. 10 and 11 into the sim- 
ulated recoil energy distribution h, in photoelectrons, are 
detailed below. 

The recoil energy spectrum obtained from the simu- 
lation is first multiplied by the £ e ff value under test to 
convert it to a spectrum with energies in keV (electron- 
equivalent) and convolved with a Gaussian energy res- 
olution with standard deviation Ry/E, where R is the 
resolution parameter under test. Next, the recoil energy 
spectrum is multiplied by the measured light yield L y 
to obtain a spectrum in photoelectrons. The number of 
photoelectrons 7V po is allowed to fluctuate according to 
a Poisson distribution. The effect of the PMT gain fluc- 
tuations is incorporated by convolving the recoil energy 
spectrum in photoelectrons with a Gaussian single pho- 
toelectron resolution with 0.6^ N pc standard deviation, 
where 0.6 is the measured mean PMT single photoelec- 
tron resolution. The measured trigger efficiency function 
discussed in section II is then applied to the recoil en- 
ergy spectrum. As mentioned earlier, the resulting recoil 
energy spectrum is divided by the simulation livetime, 
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FIG. 8: (left) Recoil energy spectrum for the 23°, 26.5°, 30°, and 34.5° scattering angles, with PSD and EJ301 energy threshold 
cuts applied, for the TOF window indicated in the TOF spectrum by the vertical dashed lines (right). The black histogram 
is the recoil energy spectrum, after subtraction of the accidental spectrum shown as an orange histogram. As a reference, the 
90% measured trigger efficiency is indicated by the vertical red dashed line. The accidental spectrum expectation is obtained 
from the TOF windows before and after the main TOF peak, as indicated in the figure by the vertical orange dashed lines. 
The accidental spectrum in the window before the peak is in agreement with the one after the peak. The purple histogram is 
the TOF spectrum where the PSD cut is chosen to select 7 ray interactions in the EJ301 neutron detectors. 
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FIG. 9: (left) Recoil energy spectrum for the 39.5°, 45°, 53°, and 120° scattering angles, with PSD and EJ301 energy threshold 
cuts applied, for the TOF window indicated in the TOF spectrum by the vertical dashed lines (right). The black histogram 
is the recoil energy spectrum, after subtraction of the accidental spectrum shown as an orange histogram. As a reference, the 
90% measured trigger efficiency is indicated by the vertical red dashed line. The accidental spectrum expectation is obtained 
from the TOF windows before and after the main TOF peak, as indicated in the figure by the vertical orange dashed lines. 
The accidental spectrum in the window before the peak is in agreement with the one after the peak. The purple histogram is 
the TOF spectrum where the PSD cut is chosen to select 7 ray interactions in the EJ301 neutron detectors. 
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FIG. 10: (left) Simulated LXe nuclear recoil energy spectrum of neutrons interacting in the active LXe volume and in the 
EJ301 neutron detector, not including accidental coincidences, for the 23°, 26.5°, 30°, and 34.5° scattering angles, and for the 
TOF window used in the analysis; and TOF spectrum (right). The solid black histogram is the total spectrum while the green 
dashed histogram is the spectrum of neutrons that elastically scatter once in the LXe active volume and maybe elsewhere. 
The red dotted histogram is the spectrum of neutrons that elastically scatter multiple times in the active volume. The blue 
histogram is the spectrum of neutrons that interact only via a single elastic scatter in the active LXe volume. 
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FIG. 11: (left) Simulated LXe nuclear recoil energy spectrum of neutrons interacting in the active LXe volume and in the 
EJ301 neutron detector, not including accidental coincidences, for the 39.5°, 45°, 53°, and 120° scattering angles, and for the 
TOF window used in the analysis; and TOF spectrum (right). The solid black histogram is the total spectrum while the green 
dashed histogram is the spectrum of neutrons that elastically scatter once in the LXe active volume and maybe elsewhere. 
The red dotted histogram is the spectrum of neutrons that elastically scatter multiple times in the active volume. The blue 
histogram is the spectrum of neutrons that interact only via a single elastic scatter in the active LXe volume. 
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FIG. 12: Fits of Monte Carlo generated recoil spectra (solid lines) to the measured recoil spectra (data points) for all measured 
scattering angles. 
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computed from the neutron generator yield at the op- 
erating conditions. Since R is left as a free parameter 
during the \ 2 minimization, any additional contribution 
to the resolution not accounted for will be absorbed in 
that parameter. 

The last step involves multiplying the simulated re- 
coil energy spectrum by an overall, energy independent 
efficiency e, taken as the same for all scattering angle 
measurements, mostly due to the EJ301 energy threshold 
cut and to the uncertainty in the neutron generator yield. 
This efficiency is computed during the x 2 minimization 
as an additional parameter for the measurement at the 
recoil energies of 8.4, 10.7 and 14.8 keV. The best fit value 
from these three measurements is taken as the efficiency 
for all measurements, while its uncertainty is taken as 
the maximum deviation in the measurements. The value 
obtained is e = 0.41^;^. The 10% relative uncertainty 
on the efficiency reflects the fact that uncertainties in the 
neutron yield from the generator are at this level. 

IV. RESULTS 




'l 10 100 

Nuclear Recoil Energy [keV] 



FIG. 13: Measured C e s values as function of nuclear recoil 
energy, together with measurements from other groups [11, 
12, 16-19]. A possible additional systematic uncertainty from 
the trigger efficiency roll-off alone, as described in the text, is 
indicated by the additional error markers. 



The £ e H values obtained for all scattering angles mea- 
sured are listed in Table I. Fig. 13 shows the results 
along with those of prior measurements at low ener- 
gies [11, 12, 16-19]. 

TABLE I: Values of C e a obtained at the 8 angles used in this 
study, together with their errors as discussed in the text. 



E nl (keV) £ eff 

23° 3.0 ± 0.6 0.088±o.oi5 

26.5° 3.9 ± 0.7 0.095±g;gi| 

30° 5.0 ± 0.8 0.098tg ^5 

34.5° 6.5 ±1.0 0.121 ±0.010 

39.5° 8.4 ± 1.3 0.139 ±0.011 

45° 10.7 ± 1.6 0.143 ±0.010 

53° 14.8 ± 1.3 0.144 ±0.009 

120° 55.2 ±8.8 0.268 ±0.013 



The total uncertainty on C c s is given by a combination 
of statistical and systematic factors with the statistical 
uncertainty taken from the fit to the data. The system- 
atic uncertainties include contributions from the spread 
in nuclear recoil energies, a\E nr , and uncertainties associ- 
ated with the 57 Co light yield, ctl , the efficiency of the 
liquid scintillator threshold cut, er e , and the positions of 
the neutron generator, o> , and of the EJ301 detectors, 
a Ts . Explicitly 



+ (^f) 2 -e 2 +(^f) 2 < + (^f) 2 <- (8) 

The change in £ e ff with nuclear recoil energy, <9£ ff / dE nl , 
is computed in closed form from a logarithmic fit to 



the measured £ e ff values versus nuclear recoil energy. 
The change in the inferred £ e g values due to the un- 
certainties in the efficiency, and the positions of the neu- 
tron generator and of the EJ301 detectors, A£ c ff/Ae, 
A£ e ff/Ar g , and A£ e ff/Ar s , respectively, were all calcu- 
lated through a discrete approximation by performing 
additional Geant4 simulations where each parameter is 
varied by small amounts. 

The largest contribution to the total uncertainty in 
£ e ff comes from the spread in recoil energy. At energies 
below 6.5 keV, this is followed by the uncertainty in the 
efficiency e. This is expected as in this region, C e g varies 
most with nuclear recoil energy, and since the mean of the 
pure single elastic recoil peak is in the roll-off of the trig- 
ger efficiency curve. For energies of 6.5 keV and above, 
the statistical uncertainty from the fit is the next-to lead- 
ing contribution to the total uncertainty. This is likely 
caused by the smaller statistics acquired for the higher 
recoil energy datasets. 

The systematic uncertainty in C e ff coming from the 
uncertainty in the trigger efficiency roll-off has been in- 
vestigated by varying the trigger efficiency function. If 
one assumes the measured trigger efficiency as the true 
efficiency, then its statistical uncertainty has a neglige- 
able effect on the inferred C c g values. However, if one 
assumes that a systematic effect is responsible for the 
discrepancy between the mesured and simulated trigger 
efficiencies and takes the simulated efficiency as the true 
efficiency, then the effect on the C c s values below 6.5 keV 
is substantial. The effect is indicated by the additional 
error indicators in in Fig. 13. 
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V. DISCUSSION 

Our results suggest that £ e ff slowly decreases with 
decreasing energy, from 0.144 ± 0.009 at 15 kcV to 
0.088+°;°" at 3 keV. The agreement at 8.4 and 10.7 keV 
with the points at 8 and 10 keV of Aprilc et. al [18] is 
excellent. Considering that the two measurements were 
performed with different LXe detectors, using different 
incident neutron energies, and at different neutron fluxes, 
reinforces the accuracy of the new measurement. Our re- 
sults below 10 keV are consistent with those of Manzur et 
al. [19], within errors, but the decreasing trend observed 
is not as pronounced. However, they are incompati- 
ble with the indirect C e g measurement of Lebedenko et. 
al [12] or with considerably smaller C c s values such as 
those suggested in [25]. 

Wc have not performed any measurement of C e g below 
3 keV to avoid the low trigger efficiency region. If C c ff 
were to decrease below 3 keV with the same logarithmic 
slope as observed between 3 and 10.7 keV, the single 
clastic recoil peak for a 2 kcV measurement, for example, 
would lie at a trigger efficiency of 45%. At this stage, the 
fraction of the neutron scattering rate which produces 



measurable signals becomes comparable to the relative 
uncertainty on the neutron generator yield. Considering 
that this LXe detector has the highest light detection 
efficiencies achieved in a LXe detector, measuring £ e ff m 
the near future at lower energies is probably impractical 
and will be subject to a considerably higher systematic 
uncertainty from the trigger efficiency roll-off. 

Given that most LXe dark matter detectors use both 
the scintillation and ionization signals to distinguish nu- 
clear recoils from backgrounds and would thus benefit 
from an improved energy calibration using both signals, 
the next step will be to simultaneously measure the scin- 
tillation and ionization yield of nuclear recoils. 
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